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We analyze the screening of static diquark sources in 2-flavor QCD and compare results with 
the screening of static quark-antiquark pairs. We show that a two quark system in a fixed color 
representation is screened at short distances like a single quark source in the same color represen- 
tation whereas at large distances the two quarks are screened independently. At high temperatures 
we observe that the relative strength of the interaction in diquark and quark-antiquark systems, 
respectively, obeys Casimir scaling. We use this result to examine the possible existence of heavy 
quark-quark bound states in the high temperature phase of QCD. We find support for the existence 
of bb states up to about 2T C while cc states are unlikely to be formed above T c . 
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I. INTRODUCTION 

Correlation functions of static quark and antiquark 
sources, represented by the Polyakov loop and its her- 
mitean conjugate, have proven to be very useful observ- 
ablcs that allow to characterize thermal modifications 
of the interaction among partons at high temperature. 
At short distances Polyakov loop correlation functions 
are sensitive to the interaction between quarks and an- 
tiquarks which is transmitted through gluon exchange. 
The thermal modification of gluons propagating in a hot 
medium is reflected in the screening of the interaction 
among static quark and antiquark sources. This allows, 
for instance, the study of the running of the QCD cou- 
pling [l[ and its temperature dependence 0, Q . 

The large distance behavior of static quark correlation 
functions is sensitive to the confining properties of the 
QCD vacuum, the breaking of the string due to pair cre- 
ation in the vacuum as well as the transition to a decon- 
fined phase at high temperature and density. While at 
low temperature the correlation functions calculated in a 
SU(3) gauge theory approach zero at large distances and 
thus reflect strict confinement of quarks, they approach 
constant values in QCD with light dynamical quarks in- 
dicating the onset of string breaking at some charac- 
teristic distance. At zero temperature string breaking 
arises through the creation of a quark-antiquark pair 
from the QCD vacuum. With increasing temperature 
more complex mechanisms involving quarks and gluons 
from the thermal medium also start playing a role. In 
the high temperature phase static quark correlation func- 
tions approach constant values in QCD as well as in pure 
gauge theories suggesting that screening of static sources 
through a gluon cloud is the dominant mechanism for 
the neutralization of the color charge of static external 
sources. 

Information deduced from the temperature depen- 
dence of the change of free energy induced in a thermal 
medium due to the presence of external sources has been 
used to model the thermal modification of the potential 
between heavy quarks. Through an analysis of the non- 



relativistic Schrodinger equation, which may be a good 
approximation for heavy quarks, it has been concluded 
that some heavy quark bound states may also survive 
the transition to the hi gh t emp erature phase of QCD 
0, H !, EI El El El 13 EI- Recently it has been 
argued that also more exotic, colored bound states made 
up of quark antiquark pairs or quark-quark pairs may 
exist in a thermal medium at high temperature [ll[. To 
quantify such a scenario for quark-quark systems requires 
a better understanding of the interaction among quarks 
at high temperature and the screening of colored diquark 
states. 

First studies of quark-quark interactions at finite tem- 
perature have been performed in quenched QCD fl6j and 
some results for QCD with light staggered fermions have 
been reported at the Lattice conference in 2005 [3 EH • 
In this paper we will focus on an analysis of static quark- 
quark correlation functions in 2-flavor QCD. We will cal- 
culate free energies of diquark systems at finite tem- 
perature and compare with results obtained for quark- 
antiquark systems at the same temperature. We will also 
calculate the induced quark number that arises from the 
presence of a diquark source that has non- vanishing trial- 
ity in a thermal bath of quarks and gluons. This provides 
information on the response of the thermal medium to 
the presence of external sources and allows to distinguish 
between the quark dominated screening mechanisms at 
low temperature (string breaking) and the gluon domi- 
nated screening at high temperature. 

This paper is organized as follows. In the next Section 
we introduce the framework for our study of diquark cor- 
relation functions, discuss the basic setup for our numeri- 
cal calculations and present some basic results on thermal 
properties of diquark and quark-antiquark free energies. 
In Section III we present results on free energies of color 
anti-triplet and sextet diquark configurations calculated 
in Coulomb gauge and discuss the gradual change from 
string breaking at low temperature to screening at high 
temperature which becomes transparent in the behavior 
of the net quark number induced in a thermal medium 
due to the presence of static quark sources. In Section IV 
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we discuss consequences of the observed Casimir scaling 
in diquark and quark-antiquark systems for the possible 
existence of colored quark-quark bound states at high 
temperature. We finally give our conclusions in Section 
V. 



II. DIQUARK FREE ENERGIES 

The general framework for the analysis of thermal 
properties of QCD in the presence of n static quarks 
and m static antiquarks has been formulated by McLer- 
ran and Svetitsky We concentrate here on 2 quark 
systems for which the lattice regularized QCD parti- 
tion function with n/ light dynamical quark degrees of 
freedom and static sources represented by the operator 

L oQ( r ) is § iven b y> 

Z% Q ((3,m,r) = J J\dU WiV L$> Q (r) (det D(m, M )f ' /4 x 

X.V 

e -PS G (U) (1) 

where U XjV € SU(3) are gauge field variables defined 
on the links of a 4-dimensional lattice of size N^N T ; 
(3 = 6/<? 2 denotes the gauge coupling, m is the bare quark 
mass for the n/ degenerate quark flavors. Furthermore, 
D(rh, /i) denotes the fermion matrix, for which we use the 
staggered fermion discretization scheme. Although all 
our calculations have been performed at vanishing quark 
chemical potential (ix) , we have indicated in Eq. [T] the 
dependence of the fermion determinant on p.. We will 
make use of this in our discussion of the quark number 
induced in a thermal medium due to the presence of ex- 
ternal sources. 

In Eq. [1] Lqq (r) denotes the operator for two static 
quark sources separated by a distance r. The superscript 
(c) indicates the color representation of the diquark sys- 
tem which we will specify later. A static quark source is 
represented by the Polyakov loop, 

P(x) = (z nf (J3,m)f T Y[ U (xchS)fi , (2) 

x = l 

and the source for antiquarks is given by the hermitean 
conjugate, pt(x). In the definition of the Polyakov 
loop we have introduced a renormalization constant, 
z„ f (/3, m), that removes divergent self energy contribu- 
tions to the static quark sources and insures that the 
diquark partition function normalized with the partition 
function in the absence of any sources, 

Z(P,m)= fj{dU x<v {detD{m^)) n '' A e-^ u \ (3) 

X.V 

has a well defined continuum limit. The renormalization 
constants for static quark sources have been determined 
previously for nj — (quenched QCD) [2(| and 2-flavor 



QCD (n f = 2) 1 for the action and for values of the 
gauge coupling used also in this study. 
The normalized QQ partition function, 

is the correlation function for two static quark sources 
separated by a distance r. Its logarithm defines the free 
energy of the static diquark system, i.e. the change in 
free energy of a thermal medium at temperature T that 
arises form the presence of two static sources, 

Fg(r,T) = -TlnCW(r,T). (5) 

We note that the discussion of diquark partition func- 
tions and free energies given so far is completely analo- 
gous to the discussion usually given for the expectation 
value of the Polyakov loop 1 , 

X 

and the related static quark free energy, 

F Q (T) = -T In (L Q ) . (7) 

The Polyakov loop operator Lq simply has been replaced 

by the diquark operator Lqq(r). 

Similarly one may, of course, also discuss the partition 
function and free energies of quark-antiquark systems. 
To be specific we introduce here the so-called color aver- 
aged operators 2 for a diquark and quark-antiquark sys- 
tem, 

L<g(r) = TrP(0)TrP(r) , (8) 
4 a J(r) = TrPt(0)TrP(r) . (9) 

Like the Polyakov loop also the operator defining the 
diquark partition function is not invariant under global 
Z(3) transformations. The static quark and diquark 
sources represent states of non-vanishing triality. Due 
to the Z(3) symmetry of the SU(3) gauge action parti- 
tion functions with non-vanishing triality will be zero in 
quenched QCD. To insure a proper definition of diquark 
partition functions also for the pure SU (3) gauge theory 
one should introduce a symmetry breaking term which 
is removed again after taking the thermodynamic limit. 
In the presence of dynamical quarks the Z(3) symmetry 
is, however, explicitly broken and the partition functions 
with non-zero triality are well defined. The fermion sec- 
tor of the QCD partition function, Eq.[31 provides the ap- 
propriate number of quark and antiquark contributions 



1 Here and everywhere else we use the normalization of traces in 
color space Trl= 1. 

2 Note that we always use renormalized operators for static quark 
sources. 
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FIG. 1: Diquark (open symbols) and quark-antiquark (filled 
symbols) free energy for color averaged sources separated at 
a distance expressed in units of the square root string tension 
r^fa at temperatures below and above the transition temper- 
ature. 



tems; the difference in Casimir factors for one gluon- 
cxchange in quark-quark and quark-antiquark systems 
suggests that the attractive interaction in an anti-triplet 
QQ system is half as strong as in a singlet QQ system. 
The color averaged free energies receive also contribu- 
tions from octet and sextet representations of a QQ and 
QQ system, respectively. These contributions are, how- 
ever, repulsive and thus are exponentially suppressed in 
Polyakov loop correlation functions at short distances. 
An additional difference in free energies of QQ and QQ 
arises from the fact that due to the non- vanishing triality 
of a QQ state the free energy of a diquark system will 
receive an additional screening contribution which is not 
present in a QQ system. 

To quantify the differences seen in Fig.[T]and to explore 
in more detail the mechanism that leads to screening of 
color charges in a hot medium it is advantageous to study 
diquark and quark-antiquark correlation functions in def- 
inite color channels which can be defined in fixed gauges. 



such that only states with vanishing triality will finally 
contribute to the partition function. We will elaborate 
further on this in Section lill Bl 

The analysis of diquark free energies presented in this 
paper has been performed on data samples that have 
been used previously to analyze free energies of static 
quark-antiquark systems Q. These data samples have 
been generated on 16 3 x 4 lattice for 2-flavor QCD using 
an improved staggered fermion action and quark masses, 
that correspond to a light pseudo-scalar (pion) mass of 
about 770 MeV. For further details on this data sample 
and further details on the action and simulation param- 
eters we refer to [3j and references quoted therein. The 
SU(3) pure gauge data used in Section IIII CI have been 
obtained on lattices of size 32 3 x 4, where we used the 
tree-level Symanzik improved gauge action j4j. For fur- 
ther details see (2lj . 

In Fig. [1] we compare the color averaged free energy 
of a static diquark system with that of a static quark- 
antiquark system both calculated in 2-flavor QCD. The 
figure shows free energies as function of the separation 
of the two sources measured in units of the zero tem- 
perature string tension, ryfa. In both cases renormal- 
ized Polyakov loop operators, Eq. ® have been used; no 
further normalization of the relative magnitude of QQ 
and QQ free energies thus is needed nor has it been per- 
formed. Apparently the free energies of a diquark and 
a quark-antiquark system approach the same large dis- 
tance limit at all temperatures. In fact, as noted previ- 
ously 0, 0| , this asymptotic value equals twice the free 
energy of a static quark represented by the renormalized 
Polyakov loop, Eq. [7J This observation confirms our pic- 
ture of screening at large distance; the two static sources 
are screened independently and the thermal bath does 
not distinguish between quark and antiquark sources. 

At short distances the QQ free energies drop faster 
than the QQ free energies. This is expected from a 
perturbative analysis of the interaction in 2 quark sys- 



III. ANTI-TRIPLET AND SEXTET DIQUARK 
FREE ENERGIES 

In addition to the color averaged diquark source in- 
troduced in Eq. [9] we will consider here sources in dif- 
ferent color channels. This is analogous to studies of 
color singlet and octet free energies performed for static 
quark-antiquark systems. A system of two quarks can 
be either in a color anti-triplet (antisymmetric) or color 
sextet (symmetric) state, 

3® 3 = 306. (10) 

The corresponding operators defining the anti-triplet and 
sextet partition functions are given by [f| 

L QQ(T) = ^TrP(0)TrP(r) - iTrP(0)F(r) , (11) 

L QQ^ = ^TrP(0)TrP(r) + iTrP(0)P(r) . (12) 

These operators are, however, gauge dependent and the 
corresponding partition functions can only be studied in 
a fixed gauge. We will analyze them in Coulomb gauge. 
The partition function for the color averaged diquark sys- 
tem can be written as a superposition of partition func- 
tions for color anti-triplet and sextet systems, 

4 a J(/?,m,r) = l Z ®(p,m,r) + h^ Q (P,rh,r) . (13) 

It has already been observed in analogous studies of 
singlet and octet free energy in quark-antiquark systems 
that the large distance behavior of free energies does 
not depend on the relative color orientation of the two 
static sources. As shown in Fig. [2] this is also reflected 
in the large distance behavior of diquark free energies. 
The same gauge invariant asymptotic value is approached 
by QQ free energies calculated with sources in the anti- 
triplet and sextet representation. 
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FIG. 2: Static diquark free energy in the anti-triplet (open 
symbols) and sextet (filled symbols) channels below (upper 
figure) and above (lower figure) the transition temperature of 
2-flavor QCD. 



A. Screening of anti-triplet diquarks 

We want to explore here in more detail the relation 
between the r-dependence of static diquark and quark- 
antiquark free energies. In leading order perturbation 
theory the interaction in a QQ singlet state is twice 
as strong as in a QQ anti-triplet state. This also car- 
ries over to finite temperature perturbation theory where 
in leading order the Coulomb interaction is modified 
only through a Debye-screening term that reflects the 
in-medium properties of the exchanged gluons, 



F, 



(3) 



QQ,pert 



(r,T) 



2 QQ,pert 

3 r 



(r,T) + const. 



-m D (T)r 



const. (14) 



This suggest that at least at high temperatures the r- 
dependent QQ and QQ interaction part in the free ener- 
gies obeys Casimir scaling. However, the gluonic screen- 
ing cloud in a QQ and QQ system will differ as a func- 
tion of distance. While we just have deduced from the 
long distance behavior of free energies shown in Fig. [2] 
that at infinite distances the screening clouds contribute 
the same amount to the free energy of a diquark and 
quark-antiquark system, respectively, this clearly has to 



be different at short distances. A QQ system in a color 
singlet state does not require a screening cloud at short 
distances; it is color neutral by construction. However, 
an anti-triplet QQ state needs to be screened even when 
the separation between both sources vanishes. In fact, 
we expect that at short distances the free energy of the 
screening cloud for an anti-triplet QQ state is identical 
to that of a single static antiquark, Fq(T). 

The above considerations suggest a simple relation be- 
tween free energies of anti-triplet QQ states and color 
singlet QQ systems, which should be valid at least at 
high temperature as well as at large and short distances, 
respectively, 



2(F®(r,T)-F Q (T)). 



(15) 



In Fig. [3] we compare the left and right hand sides of 
Eq. [15] at various temperatures. As can be seen the 
equality indeed holds very well at all distances for tem- 
peratures above the transition temperature, while some 
differences show up at intermediate distances below T c . 
Nonetheless, at short distances the equality in Eq. [15] 
holds at all temperatures which shows that indeed a well 
localized anti-triplet diquark system is screened in the 
same way as a single antiquark. Moreover, we note here 
that Eq. [15] is found to hold also in the deconfined phase 
of SU(3) pure gauge theory [18| . 



B. Screening and string breaking 

We want to elaborate here a bit more in detail on the 
differences seen in the screening of color singlet quark- 
antiquark and anti-triplet diquark free energies below T c . 
As noted before the contribution of states with at least 
one additional quark is needed in the diquark partition 
function to project on the sector of vanishing triality. 
At low temperature it is, in fact, a single quark which 
dominates the additional contributions in a diquark par- 
tition function for small separations of the two quarks. 
At large distances, on the other hand, two antiquarks 
will be needed to screen the color charges of the two well 
separated quarks. The situation in a quark antiquark sys- 
tem is quite different. No additional quarks are needed 
for screening of the sources at short distances, while a 
quark-antiquark pair is needed for screening at large dis- 
tances, i.e. the net quark number will always vanish. 

At low temperature the screening mechanism at large 
distances generally is referred to as string breaking. At 
high temperature such a screening mechanism based on 
the presence of dynamical quarks in the thermal medium 
or the vacuum becomes less relevant; screening with glu- 
ons becomes the dominant mechanism. To explore the 
features of screening and string breaking and the change 
of screening mechanisms at low and high temperature 
on a more quantitative level we have analyzed the net 
quark number induced in a thermal medium due to the 
presence of a diquark system separated by a distance r. 



5 



1.5 



1.0 



0.5 




Fi S (r,T)to" 



□ 2(F^ Q (r.T)-F Q {T))fa 1 ' 



0.87 ■ 
</2 0.90 ' — •— 

rg 0.96 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 



-0.5 



mt * m • 




T/T c 
1.002 1 — ■— 
1.02 i— •- 
1.07 -^»r- 
1.11 1 — ▼ — 
1.16 i— ♦- 
1.23 ' — • 
1.36 >— e- 
1.65 — e- 
1.98 >— »- 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 



■ Fig(r,T)fo" 
□ 2(Fg (r,T)-F Q (T))/o 1 ' 



0.5 



-0.5 



QQ 



T/T c = 0.87 
T/T = 0.90 
T/T„ = 0.96 



U ******** * * * 



0.15 
0.1 
0.05 


-0.05 
-0.1 

-0.15 
-0.2 

-0.25 



QQ 



T/T = 1.002 
T/T = 1.02 
T/T c = 1.07 
T/T c = 1.16 
T/T = 1.50 



FIG. 3: Comparison of anti-triplet diquark free energies (open 
symbols) and singlet free energies of a quark-antiquark system 
(filled symbols) below (upper figure) and above (lower figure) 
the transition temperature of 2- flavor QCD. 



FIG. 4: The net quark number induced by color averaged 
sources separated by a distance r^fo at temperatures below 
(upper plot) and above (lower plot) the transition tempera- 
ture in 2-flavor QCD. 



The net quark number can be obtained from the static 
diquark free energy, Eq. [51 by considering its parametric 
dependence on the quark chemical potential through the 
fcrmion determinant as indicated in Eq. [3J A derivative 
of In Zqq (r, T) with respect to n/T evaluated at n/T = 
then yields the net quark number, 



N$hr,T) = (N q ) t 



N q L 



&(r,7) 



(16) 



where N q is the quark number operator in 2-flavor QCD, 



N„ = -Tr 
q 2 



D _1 (r?i,0) 



dD(rh, /j.) 



(17) 



Similarly we obtain the net quark number induced by a 
single static quark source, 



Nq(T) = (N q ) Q = 



iV g TrP(0) 
TrP(0) 



(18) 



For a quark-antiquark system it is straightforward to 
see that the induced net quark number indeed is strictly 



zero for all temperatures, i.e. N^ n (r,T) = 0. However, 

WW 

as argued above, for a single quark or diquark system a 
net quark number has to be induced to generate a state 
with vanishing triality. In Fig. [¥] we show the net quark 
number induced by color averaged diquark sources sep- 
arated by a distance r. We observe a clear separation 
between the short and long distance behavior of this ob- 
servable. While for small separation of the sources a net 
quark number (Nqq > 0) is induced, a net antiquark 
number (Nqq < 0) is induced at large separations. This 
effect becomes more pronounced with decreasing temper- 
ature. Above the transition temperature, however, the 
net quark number rapidly goes to zero. 

The properties of Nqq become more apparent when we 
consider also the net quark numbers induced by static 
sources in an anti-triplet state. In Fig. [5] we show the 
values for the induced net quark number in color aver- 
aged and anti-triplet diquark systems in the large dis- 
tance limit as well as for the shortest distance avail- 
able, rT = 1/N T . At all temperatures we find that the 
net quark number induced for small separations of the 
sources in the anti-triplet channel is positive and equals 
the net quark number needed to 'neutralize' a static anti- 
quark source. At large separations, on the other hand, it 
is minus twice that needed to 'neutralize' static antiquark 
source. At short distances a QQ anti-triplet diquark 
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a static antiquark sources, Nq (upper triangle) , which equals 
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distance roa = 1/4 (squares). 
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FIG. 6: Comparison of anti-triplet and sextet diquark free en- 
ergies at short distances with the free energy of static fermions 
in triplet and sextet representations. The particular combi- 
nation shown eliminates the short distance Coulomb terms 
contributing to the diquark free energies. The sextet free en- 
ergies of static quarks have been computed in Ref. [2ll |. 



thus is 'neutralized' by a light quark from the thermal 
medium, while at large distances both quark sources are 
individually 'neutralized' by a light antiquark. While at 
finite temperature Nqq also receives contributions from 
other (multi-quark) states, this mechanism becomes in- 
creasingly dominant at low temperature. In fact, the 
data shown in Figs. [4] and [5] support the expected limit- 
ing behavior, 



lim N QQ (r,T) = 




(19) 



with r c being a typical hadronic scale, r c ~ 1.5 j \fa ~ 
0.8 fm, that characterizes the string breaking radius at 
zero temperature. 

In all cases considered here the induced net quark num- 
ber rapidly becomes small at high temperature. In fact, 
already for T>1.2 the deviations from zero are smaller 
than 0.02 at all distances. Thermal fluctuations of the 
quark number and the largely available number of ther- 
mal gluons thus seem to be sufficient to screen the static 
color sources at high temperature. 



While the large distance behavior of the sextet and anti- 
triplet free energies is identical (see Fig. [5]), we again 
expect that at short distances the sextet diquark state 
is screened differently from the anti-triplet. In fact, we 
expect that screening of a sextet diquark at short dis- 
tances is identical to screening of a single fermion in a 
color sextet representation, i.e. when combining sextet 
and anti-triplet free energies in proportion to the cor- 
responding Casimir factors such that the repulsive and 
attractive Coulombic contributions cancel we expect to 
find 



hm(2F^(r,T)+F^(r,T)) = 2F™(T) + F Q (T) (21) 

where Fq (T) is the free energy of a fermion in a 
color sextet representation given by the corresponding 
Polyakov loop expectation value, 



,(3) 



(6), 



>/T 



(22) 



where 



C. Screening of sextet diquarks 

The considerations of the previous section can also be 
carried over to sextet diquark systems. The perturbative 
one-gluon exchange among quarks in a sextet represen- 
tation is repulsive. Also this interaction is screened at 
finite temperature. At high temperature we thus expect 
to find for the genuine 'interaction part' in the sextet free 
energy, 



(r,T) 



■-F (1) -(r T) 
4 QQ' ' > 



1 a i T ) r -m D (T)r 

3 r 



(20) 



L 6 (f) = ^(TrP(f)) 2 -iTrPt (af) ( 23 ) 

is obtained from standard group theoretical relations. 
Eq. [21] is indeed fulfilled at all temperatures in 2-flavor 
QCD as well as in the deconfinement phase in SU(3) pure 
gauge theory, as can be seen in Fig. [Sj We note that the 
expectation value of a Polyakov loop in the sextet rep- 
resentation is, of course, gauge invariant. The results 
on the short distance behavior of sextet as well as anti- 
triplet quark-quark free energies thus reflect gauge in- 
variant screening properties of colored diquark systems. 
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IV. HEAVY DIQUARK BOUND STATES IN 
THE QUARK GLUON PLASMA? 

The analysis of diquark and quark- antiquark free en- 
ergies presented in the previous section suggests that 
aside from the different screening mechanisms for QQ 
and QQ systems the r-dependence of the free energies 
obeys Casimir scaling above T c ; i.e. the interaction in a 
diquark system is only half as strong as in a QQ system. 
An immediate consequence of this is that the tempera- 
ture dependent coupling defined in terms of systems, 
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bb (T) 
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TABLE I: Dissociation temperatures, Tdi s , for heavy quark- 
antiquark and diquark bound states in units of the transition 
temperature T c and dissociation energies for these states at 



3r 2 dF®(r,T) 
2 dr 



coincides with that extracted previously from color sin- 
glet free energies 

The close relation between color singlet quark- 
antiquark free energies and anti-triplet diquark free ener- 
gies, Eq.[15l of course also carries over to the correspond- 
ing energies (U) and entropies (5) deduced from the free 
energy, F = U — TS, through appropriate partial deriva- 
tives with respect to temperature at fixed separation, r, 
of the sources. Free energies as well as combinations of 
energy and entropy of static quark systems have been 
used to construct temperature dependent, effective po- 
tentials for heavy quark systems that allow to model the 
temperature dependence of heavy quark bound states in 
the high temperature phase of QCD. Although such con- 
structions are not rigorous and will include model as- 
sumptions, it seems to be clear from the analysis pre- 
sented in the previous section that any potential con- 
structed for color singlet quark-antiquark systems also 
puts stringent constraints on the interaction in a diquark 
system. We expect to have 



Here we have used parameters for the coupling in the 
(24) Coulomb term and the linear rising confinement term 
of the zero temperature potential that describe well the 
lattice results for the static quark potential, (a, a) = 
(0.385, 0.224GeV~ 2 ) and adjusted the screening mass, 
/x(T), such that the finite temperature potential re- 
produces the melting temperature of the charmonium 
ground state, J/ip, at a temperature consistent with 
lattice studies of spectral functions, Tdis(cc) ~ 1.5T C 
[22l I23I [24| . With this ansatz we find dissociation tem- 
peratures for ground states in the bottomonium system 
(T) and the corresponding diquark states given in Ta- 
ble |TJ As can be seen this analysis suggests that cc states 
are unlikely to exist in the high temperature phase of 
QCD. They dissolve already at T c . However, bb states 
dissolve only at temperatures of about 2T C . 



V QQ (r,T) = -V QQ (r,T) 



If a tightly bound QQ systems can exist at high tem- 
perature it thus is conceivable that also QQ states may 
persist in some temperature regime above T c . As we have 
seen in the previous section the color charge of such a QQ 
system will in addition be screened by a gluon cloud. 

Using Eq. [^5] we have solved the Schrodinger equation 
for diquark states formed from charm and bottom quarks 
respectively, 



2m a + —V 2 + V x (r,T) 



= M?(T) §1 



= c, b. 
(26) 

As we do not want to go into a detailed modeling of 
heavy quark bound states but rather want to discuss 
some generic features expected to hold for QQ-systems 
given some knowledge on the behavior of QQ-systems, 
we have used a generic ansatz for the heavy quark po- 
tential, which has been used previously also in studies of 
charmonium and bottomonium systems, 



V. CONCLUSIONS 

We have presented results on free energies of static di- 
quark systems. A detailed comparison with correspond- 
^25) ing results for free energies of static quark-antiquark 
systems shows that the relative magnitude of the r- 
dependent part of these free energies obeys Casimir scal- 
ing although at short distances the screening of diquark 
states is quite different from that of quark-antiquark 
states; colored diquark systems at short distances are 
screened like static quark sources in the same color rep- 
resentation. Using these results in a model for heavy 
diquark bound states we showed that it is likely that in 
a quark gluon plasma bb states can form at temperatures 
well above the transition temperature while cc states are 
unlikely to exist above T c . 



V QQ (r,T) = Ul-e-^>)-%-W 
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